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Abstract: We have analyzed the cascade of secondary electrons in diamond and 
amorphous carbon generated by the thermalisation of a single Auger electron. The 
elastic electron mean free path was calculated as a function of impact energy in the 
muffin-tin potential approximation. The inelastic scattering cross section and the en- 
ergy loss of the electron (expressed in terms of differential inverse mean free path) were 
estimated from two "optical" models, that utilise the measured dielectric constants of 
the materials. Using these data, a Monte-Carlo model describing the time evolution of 
the cascade was constructed. The results show that at most around 20 — 40 secondary 
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cascade electrons are released by a single Auger electron in a macroscopic sample of 
diamond or amorphous carbon. Consideration of the real band structure of diamond 
reduces this number further. The release of the cascade electrons happens within the 
first 100 femtoseconds after the emission of the primary Auger electron. The results 
have implications to planned experiments with femtosecond X-ray sources. 
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1 Introduction 



Radiation damage prevents structure determination of single biomolecules and other 
non-repetitive structures at high resolutions in standard electron or X-ray scattering 
experiments Cooling can slow down sample deterioration, but it cannot eliminate 
damage- induced sample movement within the time needed to complete conventional 
measurements [ffl-H]. Emerging new X-ray sources, like free-electron lasers (FEL) PLBfl, 
will offer new possibilities in imaging. Analysis of the dynamics of damage formation on 
a sample in an X-ray FEL beam suggests that the conventional damage barrier (about 
200 X-ray photons/A 2 at 12 keV energy) |2j, may be extended substantially at very 
high dose rates and very short exposure times A new dynamic barrier of radiation 

tolerance has been identified at extreme dose rates and ultra short exposure times ||- 
[7|. This barrier is several orders of magnitude higher than previous theoretical limits 
in conventional experiments. The calculations show that at these extremes, sections of 
molecular transforms from single macromolecules may be recorded without the need 
to amplify scattered radiation through Bragg reflections 

At 1 A wavelength, about nine-tenth of the interacting photons will deposit en- 
ergy into a biological sample, causing damage mainly through the photoelectric effect. 
The departing photoelectron leaves a hole in a low lying orbital, and an upper shell 
electron falls into it. This electron may either emit an X-ray photon to produce X- 
ray fluorescence or may give up its energy to another electron, which is then ejected 
from the ion as an Auger electron. The probability of fluorescence emission or Auger 
emission depends on the binding energy of the electron. In biologically relevant light 
elements, the predominant relaxation process (> 99%) is through Auger emission, and 
most photoelectric events ultimately remove two electrons from these elements (C, N, 
O, S). The two electrons have different energies and leave the atom at different times 
(for a more detailed description, see ||]8|,||). 

In very small samples (like atoms and single molecules), the primary photoelec- 
trons and the Auger electrons may escape from the sample without further interactions. 
However, in larger samples, these electrons will become trapped and thermalised. Ther- 
malisation involves inelastic electron-atom interactions, producing secondary cascade 
electrons. Here we analyze the specific contribution of Auger electrons to the ion- 
ization of a macroscopic sample through secondary cascade processes. We selected 
two different carbon compounds (diamond and amorphous carbon) as models for the 
calculations. 

Auger electrons and photoelectrons propagate through the medium in a different 
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Figure 1: Elastic mean free path (A e z or EMFP) of electrons in diamond or amorphous carbon plotted 
as a function of electron energy E. Solid line corresponds to the EMFP of electron in diamond, dotted 
line shows the EMFP of electron in amorphous carbon. 

manner. Their de Broglie wavelengths are XAuger ~ 0.8 A and Xphotoei ~ 0.1 A, respec- 
tively, and XAuger is comparable with atomic size. This implies that Auger electrons 
interact multiply with neighbouring atoms, while moving through the system of atoms 
in the solid fllQfl . Moreover, since the energy of Auger electrons is low (around 0.25 keV), 



the interaction potential must include a non-local exchange term which makes accu- 
rate description of the interaction complicated. In contrast, photoelectrons propagate 
almost freely through the medium, and their interaction with (single) atoms in the 
medium is well described by the Born approximation Therefore, in samples 



of intermediate size the low energy Auger electrons are more likely to cause significant 
ionization than the higher energy photoelectrons. The energy dependence of the mean 
free path (MFP) of electrons [O, 14" ] in carbon implies that the MFP of a photoelec- 



tron is of the order of hundred Angstroms whereas the MFP for the Auger electron 
is only a few Angstroms (> 4 A). This implies that in samples of intermediate size a 
photoelectron scatters only a few times before leaving the interaction region, while the 
Auger electron will have multiple interactions. 

In section 2 we quantify the elastic and inelastic interactions of Auger electrons 
with atoms within a solid. Using a Monte-Carlo (MC) simulation we then model the 
secondary electron cascade caused by inelastic interactions of the primary electron and 
subsequent secondary electrons with atoms. In section 3, the results of 500 computer 
simulations of different cascades are presented. These results give the estimated average 
ionization rate as a function of time. Finally, in section 4 we list our conclusions. 
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Figure 2: Energy loss function, Im[— e(q = 0,E) x ] for diamond plotted as a function of photon 
impact energy E. 

2 Secondary electron cascade in a solid 

Our study of the secondary electron effects was performed for two forms of carbon: 
diamond (p = 3.51 g/cm 3 ) and graphite-like amorphous carbon (p = 2.21 g/cm 3 ). 

Low energy electrons (E « 250 eV) may undergo elastic and inelastic collisions 
with atoms (electrons and nuclei) in a solid. Since the corresponding electron wave- 
length is comparable with atomic dimensions, multiple scattering of the electron |15| 



on neighbouring atoms have to be calculated quantum-mechanically (QM). The QM 
exchange terms must then be incorporated into the interaction potential. 



2.1 Elastic scattering 



Calculation of elastic scattering amplitudes and angular distributions can be done ac- 



curately by the partial wave expansion technique [11, 12]. In particular, the differential 
elastic cross-section —^-(E) for scattering of an electron on the atom is expressed, 



de 



using the phase shift Si of each partial wave as follows fTT|,|T2" 

da e i 



de 



c\ oo 

<E) = ^152(21 + 1)^)^(008(9)) 

K 1=0 



where k is the wave number, corresponding to the electron impact energy E, Pi(cos(6)) 
denotes the Legendre polynomial of order /, 9 is the scattering angle and the sum goes 
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Figure 3: Energy loss function, Im[— e(q = 0, E) J ] for amorphous carbon plotted as a function of 
photon impact energy E. 

over all partial wave contributions 1 = 0,... oo. The total elastic cross-section a e i(E) 
may be obtained from ([I]) by integration over 6. The corresponding elastic mean free 



path, (EMFP) X e i jl3],[T6| can be calculated as : 

X d \E) = Na el (E), 



(2) 



where N denotes the atomic density in the solid. In order to obtain the phase shifts, Si 
in ([I]), one should solve the respective radial wave equations for each partial wave with 
the approximate form of the exchange potential. To perform these calculations we used 
programs from the Barbieri/Van Hove Phase Shift package |TjJ. First we determined 
the radial charge density for a free atom, then calculated the radial muffin-tin potential 
|10|,[I5| for atoms embedded in a solid (using various approximations to the exchange 
potential), and finally derived phase shifts from the muffin-tin potential. Multiple elas- 
tic scattering within a finite cluster, provided the resulting amplitude was large enough, 
was included in the calculations. Figure |l] shows the resulting EMFP for diamond and 
amorphous carbon. For large energies, the EMFPs for both diamond and amorphous 
carbon increase linearly with electron impact energy. They also scale properly with the 

^diamond Pcarbon 



medium density, 



(cf. f[4]) 18|). With decreasing energy the 



Ael,carbon Pdiamond 

EMFPs decrease monotonically until they show oscillatory features due to interference 
between low-order scattering waves. 
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Figure 4: Inelastic mean free path (A,„ or IMFP) of electrons in diamond plotted as a function of 
electron energy E. Solid line corresponds to the IMFP calculated from Ashley's model @, dotted 
line shows the IMFP calculated from the TPP-2 model (Q. 

2.2 Inelastic scattering 



An accurate treatment of inelastic atom-electron collisions in a solid is more difficult, 
especially in the case of low energy Auger electrons when multiple scattering is impor- 
tant. In fact, a fully rigorous method for including inelastic scattering is not available 
so far. Following Fermi's work [19[], the passage of a fast charged particle was treated 
through the linear perturbation caused by its electric field in the solid. Subsequent 
developments [F2D|-[2^[] made it possible to extend the dielectric formulation in order to 
provide a more comprehensive description of quantum-mechanical effects in solids. 

Generally speaking, the linear response of a solid is described by a generalized 
dielectric constant, e(q, u), that depends both on momentum hq and frequency u. 
In quantum mechanics hu corresponds to the energy transfer of the incident charged 
particle to the solid and /iq to its momentum transfer. 



It was shown [23] that the imaginary part of the dielectric constant, Jm[— e(q, c^) -1 ], 
determines the energy loss of the test charge per unit time, dE/dt, by the formula, 
dE/dt ~ / dq J duj Im[— e(q, w) -1 ]. Therefore Im[— e(q, ou)^ 1 ] is often called the energy 



loss function (ELF). It satisfies the oscillator-strength sum rule fl8 |, that relates the 
total energy loss to an effective number of free electrons per atom, Z e ff : 



J eff 



/ dEEIml-e^E/h)- 1 } 
Jo 



(3) 
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Figure 5: Inelastic mean free path (Aj n or IMFP) of electrons in amorphous carbon plotted as a 
function of electron energy E. Solid line corresponds to the IMFP calculated from Ashley's model 
(g), dotted line shows the IMFP calculated from the TPP-2 model @. 

where ftp — J (4irn a e 2 ) / m e , n a = Nap/A is the density of atoms, Na is Avogadro's 
number, p is the density of the solid, A is the atomic weight, and E is the energy loss 
of the incoming test particle. 

The energy loss of Auger electrons in a solid is dominated by the excitation of 
plasmons. At first, we expect this behavior in metals, where conduction electrons form a 
jellium-like plasma, but not in good insulators. Nevertheless, in all solids the energy loss 
is dominated by the excitation of valence electrons to the conduction band. The excited 
electron, in turn, interacts strongly with all other valence electrons. The resultant 
eigenstate is a plasma resonance. A more familiar result of similar interactions among 
atoms in a solid is the formation of optical phonons. As expected, the plasmon interacts 
strongly with the incident Auger electron. For a more quantitative explanation, let us 
examine the dielectric function e(q, oS). It shows the importance of collective modes 
for the energy loss of charged particles. If one rewrites, e = e\ + it2 then Im[— e _1 ] = 
e 2/(ei + el)- Since e 2 is small, if e± goes to at a certain frequency u> = u P , the ELF, 
Im[— e _1 ], peaks sharply at this frequency. This corresponds to excitation of plasma 
modes of frequency Up by the incoming particle. Therefore, approximating the solid 
as a gas of free electrons, models the electron energy loss well. As the width of the 
plasma resonance and its amplitude depend on the details of the plasmon coupling and 
its decay, accurate results can be expected only from detailed simulations. 

In this paper we apply the Lindhard dielectric function approach together with 



6 



u 

CD 



0.233 - 



0.232 



0.231 - 



0.23 



0.229 



A Class B 
• Class A 



80.2 80.4 80.6 



M w (GeV) 



Figure 6: Example of an electron path in a solid. 

optical-data models. The approximation proved to work well in free-electron-like ma- 
terials where the ELF Im[— e(0, registered for incoming photons shows a dominant 



peak due to well-defined volume plasmons [18. 26 



Similarly as above, the response of the medium to a passing electron of a given 
energy hu and momentum Tiq is then described by a complex Lindhard dielectric 
function [21] e(q, oS). In general e may be a tensor but it is assumed here that the 
medium is homogeneous and isotropic. In this case, e is a scalar function which depends 
only on the magnitude of frq. The probability of an energy loss Two per unit distance 
travelled by a non-relativistic electron of energy E, i.e. the differential inverse mean 



free path (DIMFP) t(E,u) |2T|, |2|, 0, |28[ , then reads : 

1 r q + dq 



t(E,u) 



nEa 



f —Im[-e(q,u)~ 
Jq- 1 



(4) 



where Oq is the Bohr radius, and : 



g± = k 1 ± V 1 - (ftu/E) 



(5) 
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for k denoting the wave number corresponding to electron impact energy E. The 
expression for q± assumes that the energy and momentum transfer for electron moving 
in the medium is the same as for a free particle in vacuum, i.e. there is no effective mass 
assumed. Integration of the DIMFP over the allowed values of u yields the inelastic 
mean free path (IMFP) through : 

\^(E)= J dur{E,u). (6) 

It follows from that the only quantity needed to evaluate t(E,u) and \i n (E) 
is the dielectric response function e(q,u). However, most existing data on dielec- 
tric response functions were obtained from photon scattering on solids, for which the 
momentum transfer is zero. The problem is how to predict the dielectric response 
function with q > 0, knowing only its optical limit (q = 0) [p7|,p^| . For that purpose 
a phenomenological optical model approach was introduced, where Im[— e(q, is 
expressed via the convolution of Im[—e(q = O,^) -1 ] with some profile function of q 
and uj. 

The two transparent optical models we apply hereafter were chosen to give a rea- 
sonable estimate of ionization rate within the accuracy required for our model. In what 
follows we will use atomic units (h — e — m — 1) if not stated explicitly. 

The optical model by Ashley f27],|28"| includes exchange between the incident elec- 



tron and the electron in the medium modeled in analogy with the structure of the 
non-relativistic M0ller cross-section : 

i r°° 

t a {E,oj) = — — / du' J Iml-e^^y 1 ] 



2-kE J 



x {F(E, id , lu) + F(E, to 1 , E + uj' — ui ) - VF(E, uj', lo)F(E, uj', E + lu' - uj)}, (7) 

where : 



F(E,u',w) = &{lu - q 2 _/2-u>' > 0) 6(o/ + q\/2 - u > 0)— — — , (8) 



1 

' Uj(uj — LU 1 ) 

and is the step function. Substituting (||) into (|7|) one obtains P7|,P5 



t a (E,uj) = J- / duj'uj' Imi-eiO^)- 1 } 



2-kE J q 



1 1 1 



yu(uj-u]') {E + uj' - lu)(E ~ uj) - uj'){E + uj' - uj){E - uj) J 

x {0 1 (E,uj',uj) + Q 2 (E,uj',uj)), (9) 

where 0i and G2 restrict the integration region over uj' and u : 

@i{E,ui',uj) = 6(0 < uj < E/2) e{0 < uj' < 2E(lu/E - 1 + y/l - uj/E)^ , (10) 

2 (E,uj',u)) = Q{E/2 <uj < 3E/4) 6 (2ui — E < u' < 2E(uj/E -l + y/l- uj/E)j . (11) 
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The Tanuma, Powell and Penn model (TPP-2) was adopted for calculating the 
DIMFP and IMFP of electrons in a solid. We have not used the TPP-2 fit for IMFP 
calculation but derived the DIMFP, and consequently IMFP, explicitly from statistical 
approximation described in fll8| . The DIMFP tt(E,uj) yields : 

t t (E,uj) = —— / du/u/Imi-efau)- 1 ] 
Lixhj Jo 

1 



^c{uj') 2 - uj' 2 + uj 2 (^c{uj') 2 - uj' 2 + uj 2 - c(o/) 
x 0(q 2 _/2 < \jc{uj') 2 - uj 12 + to 2 - c{J) < q 2 + /2), (12) 

where c(u/) = kp(uj') 2 /3, and kp(uj') is the Fermi wave number for the free-electron 
gas with plasma frequency equal to uj' : 

k F (u') = (^y as) 

The corresponding IMFP may be obtained after integrating (|Tj|) over uj, according to 
taking into account the following restrictions : 

(g_) 2 /2 < Vc 2 -uj' 2 + uj 2 -c<{q + ) 2 /2, (14) 
E — Ep < uj, (15) 

where Ep denotes the Fermi energy (see below). In particular, restriction (0) implies 
uj' < uj. The energy loss functions for diamond and amorphous carbon used in these 
calculations are plotted in Figs. |2] and [3]. In order to obtain the ELF for diamond, we 
have used optical data for diamond |2£| , |30| (E < 35 eV) and X-ray data for scattering 
of photons on carbon |3l|] (E > 49.3 eV). The ELF in the intermediate region 35 eV 
< E < 49.3 eV was fitted in order to fulfill the oscillator-strength sum rule (Q). The 
ELF for amorphous carbon was obtained from optical data |32| in the region E < 40 
eV and from X-ray data on atomic carbon |3l| (E > 72 .4 eV). As previously, the 
ELF in the intermediate region was fitted in order to fulfill the oscillator-strength sum 
rule @. Both diamond and amorphous carbon show dominant peaks in their ELF, 
corresponding to well-defined volume plasmons |26| as expected for free-electron-like 
materials. This means that the Lindhard dielectric function approximation describes 



these two solids satisfactorily (18 

Figs. | and [5] show the IMFPs of electrons in diamond and amorphous carbon, 
calculated from ([]) and (12). The IMFPs increase monotonically with impact energy, 
however the scaling with the density of the medium is not preserved explicitly. For 
low energies (E w 50 eV), the IMFP shows a characteristic rapid increase, and for the 
TPP-2 model it becomes undefined if approaching the Fermi energy Ep (cf. (|15|)). 
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It should be stressed that the first approximations used here give an upper limit 
for the total number of secondary electrons liberated by an Auger electron. We expect 
therefore, that in reality, the number of these cascade electrons will be smaller. The 
present model treats both allotropes of carbon (diamond, an insulator and amorphous 
carbon, a conductor) as free-electron-like materials, and we model their band structure 
in a free-electron gas approximation [Kj. The Fermi energy for diamond is Ep = 28.7 
eV, and for amorphous carbon Ep = 21.1 eV as obtained from the free-electron gas 
approximation. We note that the model will give more accurate results by considering 
the real band structure of the solid. The Fermi level lies then in the middle of the band 
gap at T = K for semiconductors/insulators. 

Based on these initial results, we have constructed a model, which describes the 
time evolution of the secondary electron cascade in diamond and amorphous carbon 
(cf. ||33|| , ||34|| ). The algorithm for the MC simulation is available from authors. 

The model describes the evolution of the cascade in the approximation of inde- 
pendent non-interacting electrons, neglecting long range Coulomb interactions. The 
latter assumption holds due to the emission time scales and electron energy ranges 
relevant for the simulation. We assume that on average only one elastic or inelastic 
electron-atom scattering takes place in a cluster of size X e i^ n y An electron of energy 
E (cf. Fig. ||) enters the solid and undergoes collisions with the atoms. Depending 
on the magnitude of the respective cross- sect ions, either elastic or inelastic collisions 
occur as a stochastic process (probability of collision « g ) . In elastic collisions, 
the primary electron travels through the atomic cluster of size X e i{E) and leaves after 
time At = Xel Sz) ■ For an inelastic collision the situation gets more complicated. First, 
as previously, the electron travels through the atomic cluster of size Xi n (E). After 
time At = Xi ^^} it loses part of its energy u, and transfers it to an electron of energy 
Eq in the Fermi band (Eq < Ep). Energy Eq of the electron in the band is chosen, 
according to the Fermi density of levels at T = K |I(J (with no thermal excitations 
assumed). If the total energy E + to > Ep, the secondary electron gets excited, and it 
is emitted instantaneously when the primary electron leaves the cluster. Otherwise, if 
Eq + uj < Ep, the primary electron interacts inelastically with electrons in the Fermi 
band, loosing the part of its energy u, however, no secondary emission occurs in this 
case. The process continues until the energies of all excited electrons, including the 
primary one, fall below the Fermi barrier Ep. 

For simplicity we have assumed here that there are no thermal excitations in the 
Fermi band (T = K), and this gives an upper limit of maximal ionization. If T > 0, 
then additional low-occupied energy levels above the Fermi energy become available, so 
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the effective energy barrier becomes higher, and cascading will liberate fewer electrons 
from the Fermi band. 



3 Numerical results 



MC simulations showed that the number of cascade electrons converged after five it- 
erations in both samples. A set of 500 simulations was then performed for each of the 
two samples in order to obtain a time-dependent estimate of the number of ionizations. 
In these simulations, the energy of the primary electron was fixed at E = Ep + 250 
eV. Cascading included 1 + 5 interactions (the primary impact and 5 cascade steps). 
Figure |7| shows the results. 

For diamond the average number of ionization events after the first femtosecond 
was estimated to be ~ 6 based on Ashley's model (|9]) and ~ 7 based on the TPP-2 
model (0). The number of secondary ionizations increased with time, and it saturated 
within about 40 fs with a total of 37 electrons released at the maximal ionization of 
~ 37 events (Ashley). Saturation was slower with the TPP-2 model (100 fs), and 
the total number of cascade electrons (about 18) was about the half of those ejected 
in Ashley's model. It should be stressed that in the latter case (TPP-2) the average 
number of ionizations grew slowly with time. The same scenario held also for the 
cascades in amorphous carbon. Both Ashley's and the TPP-2 models predicted 6 — 8 
ionizations after the first femtosecond. Calculations based on Ashley's model give a 
total number of around 40 cascade electrons. These electrons were released within the 
first 10 fs, after which no more ionizations occured. Calculations based on the TPP-2 
model level out at 100 fs, and the total number of electrons released in the cascade is 
only about 23. 

The IMFP at E = 250 + Ep eV calculated from (Q) (Ashley) was larger than the 
corresponding IMFP from the TPP-2 model (12) for both diamond and amorphous 
carbon. However, the most probable energy loss at this energy is less than 60 eV in 
80 % of the CclSCS clS estimated from the integrated energy loss probability density. This 
implies that the subsequent cascade is dominated by secondary electrons of energy 60 
eV and less, and at this energy, the IMFP calculated in the TPP-2 model is larger 
than the IMFP in Ashley's model. Therefore, the number of ionizations estimated in 
Ashley's model is larger for both diamond and amorphous carbon. 

We have also plotted the maximal average ionization as the function of the electron 
impact energy (cf. Fig. ||). The total number of ionizations increases linearly with 
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impact energy in the energy range between 100 + E F eV and 300 + E F eV, as expected. 

In constructing the model, we laid emphasis on formulating a reliable description 
of the Auger electron passage through a solid. Therefore we restricted ourselves to an 
estimation of the upper limit of ionizations caused by a single Auger electron. This 
approach allowed us to use first-order approximations to model electron-solid interac- 
tions. We performed our calculations in the approximation of non-interacting electrons 
in the cascade, neglecting long range Coulomb interactions. Since the maximal num- 
ber of ions in the carbonic medium caused by a single primary Auger electron is small 
(« 20 — 40) in comparison with the total number of atoms in the sample (~ 10 9 atoms 
for (100 x 100 x 100) nm 3 cube), the approximation of neutral atoms for which the 
values for the IMFP, EMFP and DIMFP were derived is supposed to work well. This 
approach is expected to be useful for any secondary electron cascade generated by 
Auger electrons released in photoelectric events. 

Moreover, for microscopic samples one may neglect the ionization rate caused by a 
photoelectron, and then approximate the total ionization rate caused by a single photo- 
electric event by the Auger-electron ionization rate. This translates into an ionization 
rate of ~ 20 — 40 secondary electrons emitted within the first 100 femtoseconds after 
the primary electron emission in diamond and amorphous carbon. 

Finally, it should be stressed that we modeled the band structure of diamond and 
amorphous carbon, using a free-electron-gas approximation. This assumption gives an 
upper estimation of the ionization rate caused by a single Auger electron. Moreover, 
the secondary electron emission was considered in case of T = K (with no thermal 
excitations in the Fermi band), and this, again, overestimates the maximal number of 
ionizations. At T > K, the effective energy barrier becomes higher, and cascading 
will excite fewer electrons from the Fermi band than in the case of T = K. 

If one considers the real band structure of the solid, then the expected total number 
of electrons ejected in the cascade decreases further as the Fermi level lies in the middle 
of the band gap at T = K (semiconductors/insulators). The effect of including the 
real band structure on the ionization dynamics is shown in Fig. |^ for diamond, for 
which the band structure is well established [j35| -|3Tll- 

4 Conclusions 

The primary photoelectrons and the Auger electrons may escape from very small sam- 
ples, however, in larger samples, these electrons may become trapped and thermalised. 
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This process leads to additional ionization and to the deposition of further energy 
into the sample. Thermalisation involves inelastic electron-atom interactions, and pro- 
duces secondary cascade electrons on a longer time scale. In this paper, we analyzed 
the specific contribution of Auger electrons to the overall ionization of a macroscopic 
sample. The results describe the evolution of Auger-electron cascades in two model 
compounds, diamond and amorphous carbon, and show that a maximum of 20 — 40 
secondary cascade electrons may be released by a single Auger electron within the 
first 100 femtoseconds following the emission of the Auger electron. A quantitative 
description of the ionization dynamics of target samples is of crucial importance to 
practically all planned experiments at X-ray free-electron lasers, ranging from imaging 
to the creation of warm dense matter. 
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Figure 7: The average ionization rate, N e , plotted as a function of time, t, for diamond 
(upper plot) and for amorphous carbon (lower plot). The energy of the primary Auger 
electron was E = 250 + Ep eV, where Ep is the Fermi energy : Ep 29 eV for 
diamond and Ep ss 21 eV for amorphous carbon. Solid lines correspond to the average 
ionization estimated from Ashley's model, dotted lines show the average ionization 
calculated from the TPP-2 model. 
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Figure 8: Maximal ionization rate, N e ^ maX) plotted as a function of the energy of the 
primary Auger electron E = E' + E F (£" = 100, 150, 200, 250, 300 eV) for diamond 
(upper plot) and for amorphous carbon (lower plot). The corresponding Fermi energies 
are : E F 29 eV for diamond and E F m 21 eV for amorphous carbon. Solid lines 
show the maximal ionization estimated from Ashley's model, dotted lines show maximal 
ionization calculated from the TPP-2 model. 
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Figure 9: The average ionization rate, N e , plotted as a function of time, t, for diamond 
after the real band structure of diamond at T = 300 K was included into the model. 
The energy of the primary Auger electron was E = 250 + Ep eV, where Ep is the Fermi 
energy : Ep m 29 eV for diamond. The energy gap at T = 300 K equals E gap = 5.46 
eV. Solid lines correspond to the average ionization estimated from Ashley's model, 
dotted lines show the average ionization calculated from the TPP-2 model. 
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